Partially ionizing the universe by decaying particles 
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We show that UV photons produced by decaying particles can partially reionize the universe and 
explain the large optical depth observed by Wilkinson Microwave Anisotropy Probe. Together with 
UV fluxes from early formed stars and quasars, it is possible that the universe is fully ionized at 
2: < 6 and partially ionized at z > 6 as observed by Sloan Digital Sky Survey for large parameter 
space of the decaying particle. This scenario will be discriminated by future observations, especially 
by the EE polarization power spectrum of cosmic microwave background radiation. 
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I. INTRODUCTION 
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It was known for many years that the universe, which 
once became almost neutral after recombination, reion- 
ized before z ^ 5 by using Gunn-Peterson (GP) test of 
high redshift quasars. It was not until recently, however, 
that direct evidences to determine the reionization epoch 
were found. 

First evidence of the reionization epoch was obtained 
by Sloan Digital Sky Survey (SDSS) Q. The spectra 
of quasars observed by SDSS team revealed that there 
existed neutral hydrogen at the redshift z ~ 6 by GP test, 
which impUes that we have started to see the completion 
of reionization, although the amount of neutral hydrogen 
does not have to be very large to explain the observed 
GP trough. 

Recent Wilkinson Microwave Anisotropy Probe 
(WMAP) observation provided second important evi- 
dence of reionization . WMAP found that the optical 
depth due to the reionization is Top = 0.17 ± 0.04, which 
may correspond to the reionization epoch z 20. 

Therefore, it is observationally required that the reion- 
ization process was very efficient at around z 20, but 
the process did not complete until 2^6. The conven- 
tional view of the reionization process is that the uni- 
verse is reionized by UV fluxes from early formed stars 
and quasars. Many calculations have shown that the uni- 
verse can be reionized at z ~ 6 — 10 Thus, it is pos- 
sible to explain the full reionization at z = 6 observed 
by SDSS but not what WMAP observed. Recent stud- 
ies 0,0 showed that the reionization epoch can be earlier 
assuming the non-standard initial mass function (IMF) 
or unrealistically large UV fluxes from galaxies. How- 
ever, the UV fluxes from stars and quasars increase as 
time and soon exceed the critical value to fully reionize 
the universe. Therefore it is expected that the universe 
achieved the full reionization much earlier than z 6 in 
this scenario. Therefore we need to consider really ex- 
treme and perhaps unrealistic assumptions about IMF, 
escaping UV fraction from galaxies and so on, to have the 
reionization process with two stages; for example, see 0. 

Alternative to along this line, one may also consider 



other sources for ionizing hydrogens and heliums. It 
could be decaying particles, but it has been considered 
that enduring photon emission from particle decay dis- 
ables the picture of GP trough confirmed by SDSS. In 
this article, however, we will show that the optical depth 
~ 0.17 can be explained naturally by decaying particles 
with rather large parameter space of mass and lifetime; 
in some cases, the lifetime could be longer than the age 
of the universe. We follow the evolution of hydrogen, he- 
lium, and electron temperature including UV sources of 
quasars, stars, and decaying particles. Here we simply 
assume the particles decay into two photons with energy 
half of the particle mass. Notice that Hansen and Haiman 
|7j tried to explain the optical depth of Top « 0.17 by a 
particle of mass ~ 200 MeV and the lifetime ~ 4 x 10^^ 
sec, which is z « 20 reported by WMAP team 0. They 
considered the decaying particle as some sterile neutrino 
which decay into electrons (and positrons) in order to 
achieve half partially reionized era before z ~ 6. Here 
we will see that photons emitted through particle decay 
keeps much smaller partially ionization before z ~ 6, and 
explains rather large optical depth. 



II. UV SOURCES FOR IONIZATION 

We assume that the particle emits two photons with 
monochromatic energy of half mass of that particle, i.e., 
E~f — The number density of 0-particle is written 

as 



10=710(O)(1+Z)^ 



(1) 



where is the lifetime of 0-particle. If > 13.6 eV, 
those photons emitted can ionize hydrogen atoms. Then 
one can write a source term for the decaying particle as 



~dt 



T4> 



(2) 



In order to calculate how many photons are emitted, the 
abundance, mass, and the lifetime of the particle should 
be fixed. We parameterize them as free parameters by 
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fl^, E^, and r^. We will see below the allowed region in 
these parameter space. Notice that we restrict ourselves 
only for < 10^ eV. This is because E^ > MeV case 
involves e^e^ pair production, and it is too complicated. 
Moreover, E^ > 0(10) MeV may destroy light elements 
such as D, T, '^He, and ^He, spoiling successful big bang 
nucleosynthesis . 

Other than decaying particles, we also include usu- 
ally considered contributions to reionizing photons from 
quasars and stars. UV photons emitted from early 
formed stars are indeed responsible for achieving full 
ionization at z ~ 6. UV photons from quasars may 
keep the ionization fraction unity, but this contribution 
is very small, and may not be needed for that purpose 
(Of course, quasars exist and there are contributions 
from thern). We follow the argument of Fukugita and 
Kawasaki 0,1^ to follow the thermal history from z > 10'^ 
including reionization epoch, calculating the ionization 
fraction of hydrogen and helium, the photon spectrum, 
and the electron temperature (see also 0, for treat- 
ment of high energy photons from decaying particle). 

Some of the difference from Ref. 0, are adoption 
of decreasing star formation rate for smaller collapsed 
objects. In Refs.0,|3 it was assumed that the collapsed 
objects which satisfy the cooling condition become star- 
forming galaxies and a constant fraction / of baryons go 
into stars. In the present paper, we introduce the star- 
forming efficiency fstar which depends on the mass of the 
objects as 
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(Af < 3 X IO^-'Mq), 
{M > 3 X IQI^Mq), 



(3) 



where M is the mass of the collapsed object. With fstar 
the mass function of galaxies obeys the form of an em- 
pirical Schecheter function at z = 0. Since we expect 
that the UV fluxes from stars reionize the universe at 
relatively low redshift (z ~ 6), we do not have to take 
extreme values for model parameters such as the UV es- 
cape fraction Fjjv and /. We take Fjjv = 0.04^ and 
/ = 1 which give the metalicity Z ss 0.01. 

In order to estimate the UV contribution froin quasars 
we use the luminosity function in Refs.fT^ and |13^. from 
which we obtain the UV emissivity ( 1/s/cm^) at 4400Aas 
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where h is the Hubble constant in units of lOOkm/s/Mpc. 
We also assume a power low spectrum for quasars. 
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FIG. 1: Ionization histories of hydrogen (HII). We plot for 
E-j = 15 eV for = lO", 10^^ 10^^ and 10^** sec in solid 
lines from the bottom to the top (at z ~ 20), while lower and 
upper dashed lines denote for E-y — 10"' eV for — 10^* and 



10 sec, respectively. 



We set the cosmological parameters following WMAP- 
only 0: Hq = 72 km/s/Mpc, Q^n = 0.29, = I ~ 
rif, = 0.047, and as = 0.9. Apparently, we are dealing 
with ACDM universe, and the decaying particle is assume 
to be a part of cold dark matter, ^ but could be hot dark 
matter with negligible contribution to the total density 
of the universe. 



III. IONIZATION HISTORY 

We follow the evolution of the fraction of ionized hy- 
drogen and helium, in addition to electron temperature. 
Figure^shows ionization history of hydrogen for i?^ = 15 
and 10^ eV for various lifetime r^. In each case, the abun- 
dance of (/)-particles is determined in such a way that the 
optical depth reads Tgp ~ 0.17, which we define by 
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J dzuT [uHii - nHii\^^], 



(6) 



where ctt is the Thomson cross section and tihii \ is the 
number density for standard recombination. We subtract 
this term in order to estimate only the effect of reioniza- 
tion. 

In cases for E^ < 30 eV, the lifetime of 0-particle is 
allowed for many orders of magnitude. For longer lifetime 
than the age of the universe, it is the value of n^/r^ 
that concerns with the amount of photon emission, so 
ionization histories of these cases are identical to that of 
T = 10^* sec shown in Fig^ Notice that the upper limit 



Particles with very small mass such as < 10^ eV can be regarded 
t We take Fjjv = 0.005 — 0.01 for longer lifetime than the age of as cold dark matter in such a way that a scalar condensate is 

the universe for E-f = 15 eV. corehently oscillating. 
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FIG. 2; Allowed region leading to Top — 0.17. Lines show 
for Ej = 15, 30, 10^ 3 x 10^ 10^ 3 x 10^ lO", 3 x lO", and 
10^ eV from the bottom to the top. We exclude those regions 
where the photon flux exceeds the cosmic X-ray background: 
typically, Ej > 10^ eV and > lO" sec. 




FIG. 3: Ionization histories of helium (Hell). We plot for 
E-j = 15 eV for = lO", 10^^ 10^^ and 10^** sec in solid 
lines from right to left (near z ~ 20), while lower and upper 
dashed hnes denote for Ej = 10^ eV for = lO" and 10^^ 
sec, respectively. 



of the lifetime is around 4 x 10^^ sec, when (/)-particles is 
(cold) dark matter, i.e., fl^ ~ flm- For shorter hfetime, 
the abundance should be much less than that contributes 
to the energy density of the universe. See Fig. [21 

On the other hand, if larger energy of photons are 
emitted, we have to be careful about cosmic X-ray back- 
grounds. We adopt the following simple formulas to com- 
pare with the calculated photon flux 11] : 
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where F^_obs is measured in units of (cm^ sr sec)^^. 

We find that long lifetime is not allowed, because the 
emitted photon flux exceeds the observed diffuse photon 
background. Thus, the lifetime is rather restricted in 
the narrow range, and should typically be < 10^^ sec. 
This constraint is taken into account in Fig. [21 thus the 
parameter space is restricted for > 10^ eV. Hydrogen 
ionizing histories for E^ — 10^ eV are shown in Fig. ^ as 
an example. They look very similar to those for E^^ — 
15 — 30 eV. In this case, even = 10^^ sec seems to be 
excluded as is shown in Fig. [21 

Evolutions of ionization fraction of hydrogen look more 
or less similar in any emitted photon energy for the same 
lifetime. In order to distinguish them, we need to look at 
the ionization history of helium. Since larger photon en- 
ergy can excite helium atoms, one can see the qualitative 
difference from low energy emitted photons. In Figs. [31 
and^ ionization histories of Hell and Helll are shown, 
respectively. We only plotted for E^ — \b and 10^ eV 
cases as in Fig. 
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FIG. 4: Ionization histories of helium (Helll). We plot for 
E^ = 15 eV for = lO", 10^^ and 10^* sec in solid 

lines from right to left (near z ~ 20), while lower and upper 
dashed hnes denote for E^ = 10^ eV for = 10^* and 10^^ 
sec, respectively. 



It is trivial to say that rather high ionization fraction 
begins at high redshift in large E^ cases. Therefore, al- 
though it will be very difficult to observe ionization frac- 
tion of helium at z > 20, it might be one of the tools for 
(dis) proving the particle decay scenario. 



IV. COSMIC MICROWAVE BACKGROUND 
SPECTRUM 

The evolutions of ionization fraction of hydrogen and 
helium look very different from those instantaneous 
reionization history usually considered in the literatures, 
although we obtain the optical depth of ^ 0.17. So it 
is interesting to know what we can see, or more pre- 
cisely, whether we can distinguish them, in the cosmic 
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FIG. 5: TE spectrum for various ionization histories. We 
plot for = 15 eV for = lO", 10^^ lO", and 10^** 
sec from the bottom to the top in thick solid lines, while 
lower and upper thick dashed lines denote for = 10^ eV 
for = 10^^ and 10^^ sec, respectively. We plot no and 
instantaneous reionization cases in thin solid lines. Higher 
optical depth case is also shown in dotted line {E-y = 15 eV, 
= 10^*^ sec, and Top ~ 0.43). In addition, WMAP data is 
placed. 



FIG. 6: EE spectrum for various ionization histories. We 
plot for E^ = 15 eV for = lO", 10^^ 10^^ and 10^* sec 
from the bottom to the top (at z ^ 10) in thick solid lines, 
while lower and upper thick dashed lines denote for E-^ = 10^ 
eV for = 10^* and 10^^ sec, respectively. We plot no and 
instantaneous reionization casea in thin solid lines. Higher 
optical depth case is also shown in dotted line {E^ — 15 eV, 
r0 = 10^^ sec, and Top ~ 0.43). 



microwave background (CMB) power spectrum. Figure 
shows the cross correlation spectra between tempera- 
ture and E-mode polarization. Our models seem to be 
less power in the small £ region than the model with 
instantaneous reionization. This is because particle de- 
cay takes place for longer period, so the power spreads to 
larger i for fixed value of the optical depth. From the TE 
spectrum, it is not apparent to distinguish between these 
models within these error bars. Even difference between 
the instantaneous reionization and ionization by decaying 
particles can hardly seen. Moreover, notice that larger 
optical depth Top — 0.43 may account for the observa- 
tional results of the TE spectrum, too. Of course, it is 
known that larger optical depth makes the first acous- 
tic peak of the TT spectrum lower. However it can be 
adjusted by introducing larger spectral index, so-called 
Top — n degeneracy. 

Unlike TE spectra, however, EE spectra seem to serve 
us a much powerful tool for discriminating among differ- 
ent ionization histories,^ as shown in Fig.|Sl Peculiar fea- 
ture of enduring photon emission from particle decays is 
larger power around £ — {a few) x 10. In addition, there is 
different feature between long and short lifetime: for long 
lifetime, the spectrum is rather flat, while there is a dip at 
£ ^ 10 for short lifetime. Thus, EE spectrum hoped to be 
seen by Planck satellite or even WMAP should (dis)prove 
the photons from decaying particles in the near future. 



§ Even in more conventional cases, differences of EE spectra seem 
to h ave abi lity to distinguish among different ionization histo- 
ries [TllTl. 



To end this section, we comment on the Compton yc- 
parameter. We also follow its evolution, and obtained 
that the final value is (several) x 10~^ in all the parame- 
ters we choose, which is well below the current observa- 
tional limit [ig- 

V. CONCLUSION AND DISCUSSION 

We have studied comprehensively the ionization his- 
tory of the universe by photons emitted from decaying 
particles, together with the effects of stars, quasars, and 
shown that the decaying particle can explain rather large 
optical depth observed by WMAP without conflicting 
other astrophysical and cosmological constraints. One of 
the amazing fact is that the lifetime of order 10^^ sec (in 
other words, z ~ 20) is not necessarily required for suc- 
cessful scenario, and larger parameter space is allowed. 
Since the early star formation ejecting UV photons may 
be difficult to lead large optical depth without adjusting 
parameters beyond their natural values, this alternative 
is at least as good as them. 

For E-y — 10^-10^ eV, larger lifetime is prohibited, be- 
cause the emitted flux exceeds the observed diffuse pho- 
ton background, and ~ IQ-'^^ — IQ-'^^ sec is only allowed 
in these cases. On the other hand, if some bump is ob- 
served in the photon spectrum, it could be a great clue 
for decaying particle scenario. 

Longer lifetime is also allowed for < 30 eV. In the 
limiting case, the decaying particle can be the dark mat- 
ter, provided that its lifetime is around 4 x 10^^ sec. No- 
tice that our situation is different from Ref. 0| , in which 
they investigated for decaying neutrino dark matter and 
thus using hot dark matter scenario. In this case, ioniza- 
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tion fraction becomes unity at z ~ 20— 30, which conflicts 
with SDSS results. In addition, EURD observation may 
falsify decaying particle scenario with ^ 15 eV, but 
it needs careful look, since their conclusion will change 
if the abundance of galactic dark matter near the Sun is 
one third of their applied value 

Observation of CMB spectrum is very powerful. Al- 
though instantaneous reionization suggests optical depth 
Top ^ 0.17, higher optical depth may be favored for the 
observed TE spectrum by WMAP if the decaying particle 
scenario is correct, where slow reionization takes place. 
It is somewhat difficult to discriminate from each other 
only by TE spectrum. However, the differences in the EE 
spectrum is tremendous so we hope that even the lifetime 
of the decaying particle will be measured by Planck. 

Finally, we comment on what can be a candidate for 
this decaying particle. In the context of no-scale super- 
gravity, scalar particles remain massless at the tree level. 
They may acquire their masses at one-loop order. Among 



them, saxion, the scalar partner of axion in supersym- 
metric theory, might be the best candidate. If we con- 
sider the hadronic axion model, the lifetime is estimated 
as [13 

where Fpq is the axion decay constant. It perfectly fits 
into the allowed region shown in Fig. |21 
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